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ABSTRACT 

The purpose of this study was to develop a model for 
the distribution of authorships — based on the initial hypothesis that 
the distribution of authorships follows a shifted Waring 
distribution — and to test the derived model and some other discrete 
probability models for goodness-of-fit against empirical data. 
Bibliographic data from 15 abstracting journals covering the 
literature in six fields — engineering, medical, physical, 
mathematical and social sciences, and humanities — were used in 
testing the goodness-of-fit of the shifted Waring distribution and 13 
other discrete probability models. The preliminary findings presented 
here are based on 60 data sets collected from 10 abstracting journals 
covering the literature in the mathematical and social sciences and 
humanities. They indicate that the promising models for the 
distribution of authorships are the shifted Waring, shifted 
generalized negative binomial, shifted negative binomial, shifted 
generalized Poisson, and shifted inverse Gaussian-Poisson 
distributions. Three advantages and possible practical applications 
of a model for the distribution of authorships include: (1) ability 
to summarize the entire frequency distribution by a few parameters of 
the model; (2) estimation of the number of entries in an author 
index; and (3) usefulness in a simulation study designed to 
determine, subject to space constraints, the maximum number of 
authors per paper to be included in an author index. Forms of the 
discrete probability models ar£ appended. (5 tables, 12 references) 
(Author/CGD) 
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ABSTRACT 

The purpose of the study is to develop a model 
for the distribution of authorships • The initial hypothesis is that 
the distribution of authorships follows a shifted Waring 
distribution. Bibliographic data from 15 abstracting journals 
covering the literature in 6 fields (engineering sciences, medical 
sciences, physical sciences, mathematical sciences, social sciences, 
and humanities) are to be used in testing the goodness-of-f it of the 
shifted Waring distribution and 13 other discrete probability 
models. The findings presented here are based on data collected from 
only 10 abstracting journals covering the literature in the 
mathematical sciences, social sciences, and humanities. The shifted 
Waring, shifted generalized negative binomial, shifted negative 
binomial, shifted generalized Poisson and shifted generalized 
Gaussian-Poisson distributions performed fairly well in the Chi 
Square tests. 
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1. INTRODUCTION 

Bibliometricians have developed models for various social 
phenomena . The three most common models are Bradford 1 s law (1934) 
for journal productivity, Lotka f s law (1926) for scientific 
productivity, and Zipf f s law (1949) for frequency of words in text. 
Models that have been established for other social phenomena are 
listed in articles by Simon (1955), Kendall (1961), Haitun (1982a), 
and Chen and Leimkuhler (1986) • 

One social phenomenon for which no model has been developed is 
collaboration in research, though a few have been suggested without 
their goodness-of-f it being tested. Just as the number of papers 
published in scholarly journals is taken as an indication of 
productivity, co-authorships of papers in journal is often taken as 
an indication of collaboration in research. 

Price and Beaver (1966) , using memos circulated among members of 
an invisible college, inferred a Poisson model for the distribution 
of authorships while Haitun (1982b) , using the same data, 
classified the distribution of authorships among those stationary 
scientometric distributions that cannot be approximated by the Zipf 
distribution . 

Worthen (1978) observed that for drug literatures the number of 
papers with one author makes up about 1/3 of the literature, those 
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with two authors about 1/3 of the remaining literature, those with 
three, 1/3 of the remainder, and so on, Goffman and Warren (1980), 
on the other hand, observed that for the schistosomiasis literature, 
the number of publications with one author makes up about 2/3 of the 
literature, those with two authors about 2/3 of the remaining 
literature, those with three, 2/3 of the remainder, and so on. The 
two observations implied a geometric model for the distribution of 
authorships . 

The objectives of this research, therefore, are to derive a 
theoretical model for the distribution of authorships and to test 
the derived model and some other discrete probability models (for 
goodness-of-fit) against empirical data • 

2. PROBABILITY MODELS 

2.1 Theoretical Model 

The initial theoretical model being proposed is the shifted 
Waring distribution. Its derivation is as follows : 

Assume that a researcher completes a project alone if he can. 
However, if he cannot, he brings in additional researchers, one at a 
time, until the project is completed. Let an attempt by the 
researcher (s) to complete a project without bringing in an 
additional researcher represent a trial. The probability that the 
project will be completed on any trial is constant. The distribution 
of the number of trials or the number of researchers required before 
the project is completed can then be described by the geometric 

2 
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distribution. However, the probability of completing the project at 
any trial without bringing in an additional researcher, i.e., p, 
would vary from project to project (types of projects include 
review, opinion on a topic, experimental research, theoretical work, 
etc) . Thus, 

P(X = r /p) = p (l-p) 3 ^ 1 ; r = 1,2,... 
For analytical convenience, we assume that p varies as a beta 
distribution with parameters a,)8, i.e., p has a density function 



f (P) = 



p*" 1 (I-P)/ 3 - 1 



B(a,jS) 



Hence , 
P(X 



B(a,jS) 

0 otherwise 
where a > 0 , p > 0 are constants and 
1 

x a-l (i_ X )0-l dx 
0 

= r(a) T( (3 ) 

r) = E [P(X = r /p)] 

P(X = r /p) f (p) dp 



if 0 < p < l 



pfl-pjr- 1 p"" 1 (l-p) ^ 1 dp 
B(a,/3) 



B(a,/3) 



p a (l-p) r+ ^- 2 dp 



5 



= 1 B(a+l,r+)S-l) 
B(a, fi) 

T(a+B) rfa+1) rfr+fl-1) 

r()S) r(a+)8+r) 

a+0 T(a+)S+r) r(0) 

JL 

(k+g -2) ; a > 0, 0 > 0, r = 1,2,... 
(a+0+k-l) 

This distribution is known as the shifted Waring, 

2.2 Empirical Models 

In addition to the shifted Waring distribution, other discrete 
distributions that will be tested are : 

1. Zipf 

2. Mandelbrot 

3 . Geometric 

4. Shifted Poisson 

5. Shifted Generalized Poisson 

6 . Logarithmic 

7 . Borel -Tanner 

8. Shifted Yule 

9. Shifted Generalized Waring 

10. Shifted Inverse Gaussian-Poisson (IGP) 

11. Shifted Binomial 

12. Shifted Negative Binomial 

13. Shifted Generalized Negative Binomial 

The form of each model and the method of estimation of its 

4 
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parameters are given in the Appendix. 

The Chi Square test will be used to determine the goodness-of-f it 
of the models at 5% and 1% levels of significance. The Chi Square 
test was chosen in preference to the Kolmogorov-Smirnov test because 
the latter is strictly applicable only if the data are continuously 
distributed (Conover, 1980) . 



3. DATA COLLECTION 

Most bibliometric studies of research collabore* ion have used 
either scholarly journals or abstracting journals or both as their 
sources of data. Because of their comprehensive nature, abstracting 
journals will be used as the source of data for this study. However, 
only research papers published in scholarly journals or presented at 
conferences will be counted. Books, dissertations, theses, and 
research papers by anonymous or corporate authors are being excluded 
because : 

(1) Some books are made up of several articles, which are at times 
written by different authors, and not all the abstracting journals 
cover monographic literature. However, book chapters or parts of 
books abstracted separately will be counted. 

(2) Dissertations and theses are of single authorships by nature. 

(3) One cannot easily ascertain the number of people responsible for 
research works credited to anonymous or corporate authors. Also, 

5 
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editors, compilers, translators, and moderators of symposia will not 
be regarded as co-authors since they are not intellectually 
responsible for the research works attributed to them. 

In order to see how well the probability models fit the 
distribution of authorships in various fields, data will be 
collected from engineering sciences, medical sciences, physical 
sciences, mathematical sciences, social sciences, and humanities. 

Engineering Index abstracts literature from all the disciplines 
in engineering sciences while Index Medicus does the same for the 
disciplines in medical sciences. However, no single abstracting 
journal covers all the disciplines in either the physical sciences, 
mathematical sciences, social sciences, or humanities. Hence, a few 
abstracting journals have been chosen to represent each of these 
fields (papers from the abstracting journals representing each field 
will not be added together because of journal overlap). 
Availability, scope, and prominence are the factors taken into 
consideration in choosing the abstracting journals. The list of the 
chosen abstracting journals and the range of the number of records 
abstracted per year by each journal are given below in Table 1. 

Six annual cumulative issues of each abstracting journal will be 
used in the study. Because some of the chosen abstracting journals 
came into being late in the fifties or early in the sixties, the 
years chosen by systematic sampling are 1961, 1966, 1971, 1976, 
1981, and 1986. However, data for philosophy will be collected for 
1967/68, 1969, 1971, 1976, 1981 and 1986 because the first issue of 
Philosophers Index was published in 1967/68. 

6 
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For each of Engineering Index, Index Medicus, Chemical Abstracts, 
Biological Abstracts, Physics Abstracts and Psychological Abstracts, 
a few issues will be selected at random for each year. For the 
others, which index relatively small number of records per year, all 
the issues will be used. 

Articles contained in each annual cumulative issue of each 
abstracting journal will be sorted into different levels of 
authorships. The number of articles at each level will then be 
counted and the resulting distribution of authorships taken as one 
data set. Hence, there will be a to^al of 90 data sets to be used in 
the analysis. 
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Table 1. List of the abstracting journals selected and the range 
of the number of records abstracted per year by each 
j ournal 



Abstracting Journal 



Field 



Number of Records 
37000 - 131000 

110000 - 400000 

21000 - 130000 

* 

145000 - 240000 
87000 - 235000 
12000 - 40000 
800 - 4000 

3000 - 4000 
2000 - 11000 
1500 - 6500 
2000 - 7000 
14000 - 35000 
3000 - 23000 
2500 - 4000 
2000 - 7000 



1. Engineering 
Index (1917- 

2. Index Medicus 
(1960- 

3 . Physics 
Abstracts (1903- 

4 • Chemical 

Abstracts (1924- 

5. Biological 
Abstracts (1926/27- 

6. Mathematical 
Reviews (1940- 

7. Statistical Theory 
& Methodology 
Abstracts (1960- 

8 . Computer 
Abstracts (1957- 

9. Sociological 
Abstracts (1953- 

10. Inter. Pol. Sc. 
Abstracts ( 1951- 

11. Economics 
Abstracts (1953- 

12. Psychological 
Abstracts (1922- 

13. Historical 
Abstracts (1955- 

14. Abstracts of 
English Studies (1958- 

15 . Philosopher 1 s 
Index (1967/68- 



Engineering 
Sciences 

Medical 
Sciences 

Physical 
Sciences 

Physical 
Sciences 

Physical 
Sciences 

Mathematical 
Sciences 

Mathematical 
Sciences 



Mathematical 
Sciences 

Social 
Sciences 

Social 
Sciences 

Social 
Sciences 

Social 
Sciences 

Humanities 



Humanities 



Humanities 
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4. RESULTS OF THE PILOT STUDY. 



At present, data have been collected from only 10 abstracting 
journals covering the literature in the mathematical sciences, 
social sciences, and humanities. Hence, the findings presented below 
are based on only 60 data sets. 

In the humanities, many of the models performed very well with 
shifted generalized Poisson and shifted inverse Gaussian-Poisson 
doing exceptionally well (see Table 2) . 

In the social sciences, only shifted generalized Poisson, shifted 
generalized inverse Gaussian-Poisson and shifted negative binomial 
distributions did very well while shifted Waring and shifted 
generalized negative binomial distributions performed fairly well 
(see Table 3) . 

In the mathematical sciences, all the models performed poorly 
(see Table 4) . 

Overall, only shifted generalized Poisson, shifted Waring, 
shifted generalized inverse Gaussian-Poisson, shifted negative 
binomial and shifted generalized negative binomial distributions 
performed fairly well in the Chi Square tests (see Table 5) . 
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Table 2 : Results of the Chi Square goodness-of-f it tests for 
the Humanities (18 data sets) 



fioaei 


no. or valid 
Results 


NO. Of 

5% 


Pass' js 
1% 


XT _ ££ T"> x. 

No. of Best 
Fits Provide* 




lb 






1 


Mandelbrot 


18 


14 


16 


4 


Geometric 


18 


6 


10 


0 


oniiuea Jroisson 


18 


7 


/ 


1 


on « ben. roisson 


18 


T *"7 
1 / 


1 / 


1 


Logariunitiic 


18 


7 


7 


0 


Dorei-i anner 


18 


1 


9 


0 


oiiiiuea xu >.e 


1 Q 

lo 


1U 


lo 


3 


Shifted Wa finer 


17 


X** 


1 7 


n 
u 


Sh. Gen* Waring 


3 


0 


0 


0 


Shifted IGP 


17 


16 


17 


6 


Shifted Binomial 


1 


1 


1 


0 


Sh* Neg. Binomial 


17 


14 


16 


0 


Sh. Gen. Neg. Bin. 


14 


13 


14 


2 



Note : 1) "No. of valid results" refers to the number of data sets 
in which the assumptions of the model are not violated. 

2) "No. of passes" refers to the number of times that there 

was no significant difference between the observed 
distribution and the expected distribution at the 
specified level of significance. 

3) "No. of best fits provided" refers to the number of times 

that the model had the highest probability level for the 
Chi Square value 
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Table 3 : Results of the Chi Square goodness-of-f it tests for 
the Social Sciences (24 data sets) 





I h U • ui vaiiu 

Results 


5% 


1% 


Fits 


Zipf 


24 


2 


2 


0 


Mandelbrot 


24 


5 


6 


0 


Geometric 


24 


5 


7 


0 


Shifted Poisson 


24 


3 


4 


1 


Sh. Gen. Pcisson 


,24 


18 


22 


2 


Logarithmic 


24 


6 


8 


0 


Borel-Tanner 


24 


6 


7 


0 


Shifted Yule 


24 


5 


5 


1 


Shifted Waring 


12 


10 


12 


2 


Sh. Gen. Waring 


13 


1 


2 


0 


Shifted IGP 


24 


19 


24 


13 


Shifted Binomial 


0 


0 


0 


0 


Sh. Neg. Binomial 


24 


17 


20 


3 


Sh. Gen. Neg. Bin. 


14 


10 


13 


2 
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Table 4 : Results of the Chi Square goodness-of-f it tests for 
the Mathematical Sciences (18 data sets) 



Model No. of Valid No. of Passes No. of Best 

Results 5% 1% Fits Provided 





18 


n 


n 


o 


Mand^l V>Tot 


18 


o 


o 


o 


vjrfctwlUfc: Ui_ lb 


1 ft 


A 
*x 


—/ 


i 






O 


A 


n 




18 


7 


g 


2 


Logarithmic 


18 


0 


0 


0 


Borel-Tanner 


18 


0 


0 


0 


Shifted Yule 


18 


0 


0 


0 


Shifted Waring 


0 


0 


0 


0 


Sh. Gen. Waring 


4 


0 


0 


0 


Shifted IGP 


14 


5 


6 


13 


Shifted Binomial 


4 


3 


3 


0 


Sh. Neg. Binomial 


14 


3 


5 


0 


Sh. Gen. Neg. Bin. 


5 


4 


5 


2 
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Table 5 : Results of the Chi Square goodness-of-f it tests for 
the three fields combined together (60 data sets) 





iiUt \J±. vaiiu 

Results 


ii u ♦ 

5% 


rabbcb 

i% 


Fits 


Zipf 


60 


11 


14 


1 


Mandelbrot 


60 


19 


22 


4 


Geometric 


60 


15 


22 


1 


Shifted Poisson 


60 


12 


15 


2 


Sh. Gen* Poisson 


60 


42 


48 


5 


Logarithmic 


60 


13 


15 


0 


Borel-Tanner 


60 


13 


16 


0 


Shifted Yule 


60 


15 


13 


4 


Shifted Waring 


29 


24 


29 


2 


Sh* Gen* Waring 


20 


1 


2 


0 


Shifted IGP 


55 


40 


47 


32 


Shifted Binomial 


5 


4 


4 


0 


Sh. Neg. Binomial 


55 


34 


41 


3 


Sh. Gen. Neg. Bin. 


33 


27 


32 


6 
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5. CONCLUSION 



Judging from the analysis of the 60 data sets collected from 10 
abstracting journals covering the literature in the mathematical 
sciences, social sciences, and humanities, the promising models for 
the distribution of authorships are the shifted generalized Poisson, 
shifted Waring (the proposed theoretical model) , shifted inverse 
Gaussian-Poisson, shifted negative binomial and shifted generalized 
negative binomial distributions. However, we have to wait for the 
results of the analysis of the remaining 30 data sets, which are to 
be collected from 5 abstracting journals covering the literature in 
the engineering sciences, medical sciences and physical sciences, 
before determining the best model for the distribution of 
authorships. 

The advantages and possible practical applications of a model 
for the distribution of authorships include : 

(i) Ability to summarize the entire frequency distribution by a few 
parameters of the model. 

(ii) Estimation of the number of entries in an author index : 
Abstracting services often add extra entries in case of multiple 
authorships. For large abstracting journals and indexes, it would be 
time consuming to determine ahead of publication how many author 
entries would be made for N, a fixed number, papers. So, in order to 
estimate the number of entries, an abstractor may randomly sample n 
out of the N papers and estimate the parameters of the model from 

14 
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the sample. The values of the parameters could then be used to 
estimate for the N papers, the number of papers with one author, two 
authors, three authors, and so on. The abstractor then should be 
able to estimate the total number of entries, and, hence, have a 
rough estimate of the size of the author index. 

(iii) Usefulness in a simulation study designed to determine, 
subject to space constraints, the maximum number of authors per 
paper to be included in an author index. Only the first author would 
be included for papers with more than the maximum number of authors 
per paper allowed. 
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APPENDIX : FORMS OF THE DISCRETE PROBABILITY MODELS 

In the following , we write 

f(x) = frequency of occurrence of x 

n = sample size 

MLE = maximum likelihood estimate 

x" = sample mean 

s 2 = sample variance 

1. Zipf 

p(x) = k x" b 

where 0 < k < 1, b >0 for x = 1,2,.... 
The maximum likelihood estimator of b satisfies the equation 
f(x) In x / S f(x) = -j'(b) / <J (b) , 

where ^ ( • ) denotes the Riemahn zeta function ; and 
k = [S x" b J" 1 

2 . Mandelbrot 

p(x) = k (x+c)- b ; 
n 

where k>0, b>0, -l<c<«> for x = 1,2,... 
The parameters are estimated by non-linear least-squares 
regression method. 

3 . Geometric 

p(x) = p (l-p)^ 1 ; 
where 0 < p < l for x = 1,2, .. . 

16 
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The MLE of p is p = 1/x 



4. Shifted Poisson 

p(x) = e" X A* -1 ; 
(x-1) ! 

where A > 0 for x = 1,2,... 

A 

MLE of A is A = x - 1 

5. Shifted Generalized Poisson 

p(x) = X,f,\, + fx-pA^*- 2 e-( X > + (x-l)A*) ; 

(x-1) ! 

where > 0, | | < 1 for x = 1,2,... 
The moment estimates of X i and X^are 

A 

V = 1 - ((x-l)/s 2 ) ^ and 
\ = (x-1) (1 -A a ) 
The maximum likelihood estimator of A^is given by the root 
of the equation 



f (x) (x-2) - n(Y-l) = 0 

and 

V A 

\ = (x-1) (1 -Aj 



6. Logarithmic Series 

p(x) = ae x ; 
x 

where 0<9< 1 , a = - [ln(l- 6)]" 1 for x = 1,2, 
The moment estimate of 6 is 

9=1- x 

(s 2 + x 2 ) 

17 



- ERIC 



.19 



The maximum likelihood estimator of e is given by the root 

of the equation 

x + 9 = 0 ; 

(1 - §) In (1 - §) 

7. Borel - Tanner 

p(x) = x x ~ 2 a*" 1 e~ ax ; 
(x-1) ! 

where a > 0 for x = 1,2,... 
The moment estimate of a is given by 
a = 1 - 1/x 

8. Shifted Yule 

p(x) = a T(x) r(a+l) ; 

r(oc+x+i) 

where a > 0 for x = 1,2,... 
The moment estimate of a is given by 
a = 1 + l/(x-l) = x/(x-l) 

9. Shifted Waring 

p(x) a rra+fl+1) T(x+B-1) ; 

a + p r(a+0+x)r()8) 

where a > 0, p > 0 for x = 1,2,... 

The moment estimates of a and /J are given by 

ot = 2s 2 ; and 

(s 2 - x(x-l)) 

j8 = (a-1) (x-1) 

Note that the mean exists only if a > 1 while the variance 

18 
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exists only if a > 2. 

The maximum likelihood estimator of a is given by the root 
of the equation 

n(x-p (a-1) - / f (k) ^ 1 = 0 ; and 

a[x(a-l)+l] [x(a-l) + x] 

P = (a-1) (x-1) 



10. Shifted Generalized Waring 

P ( x) = rfv+a) r(x+v-i) rfx+£ -p ; 

B(a,£) r(v) r(x+v+a+0-l) (x-1) ! 

where a > 0, p > 0, v > 0 for x = 1,2,... 

The three parameters can be estmated from the first three 

ascending factorial moments given by 

M[i] £v ; a > 1 

a-1 

M[2] = flv (v+1) ffl+11 ; a > 2 and 
(a-1) (a-2) 

M[3] = fly (v+lWv+2H/?+lHfl+2) ; a > 3 

(a-1) (a-2) (a-3) 

11. Shifted Inverse Gaussian-Poisson 

p(x) = (2a/?r)^ e^ 1 " 6 ) (hae) *" 1 K x _ 3/2 (a) 

(x-1) ! 

where 0<8<1, a>0, K v (z) is the modified Bessel 
function of the second kind of order v with argument z foi 
X 1 j 2 f • • • 

Estimation of parameters : 

(a) If the observed frequency distribution is 
reversed J-shape, with a large proportionate frequency in 

19 
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the unit cell. 

e 1 - P -In v(x=l) ~~\ 2 ; and 

I 2(x-l) + In p(x=l)j 

A A 1 

a = 2fx-l) (1 - 6) -5 
e 

(b) I observed frequency distribution is 

unimodal . 

The moment estimates of e and a are given by 
e = 1 - r 2s 2 - "ll -1 ; and 



1 - p 2s 2 - l -1 ; 
lJx-1) J 



a = 2(x-l) (l - 9) ^ 
9 

The maximum likelihood estimator of a is 
given by the root of the equation 

)~*][1 + (x-l)/w] - (1/n) |f f(x > R x-3/2 («> = 0 

where e = x - 1 , w = (e 2 +a 2 )^ - e , R v (z) = K v+1 (z) 

K v (z) 



and 



9 = -2e 2 ± 2e(e 2 +a 2 
a 2 



12. Shifted Binomial 

P(x) = ^ v^ p^-1 (1 _ p) v-(x-l) . 

where 0<p<l,v>0 for x = 1,2, ..,v,v+l 

The moment estimates of p and v are given by 

•t 

p = 1 - s 2 ; and 

x-l 

20 
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The maximum likelihood estimator of v is given by the root 
of the equation 

*^ '(v-j+l)" 1 fj + n In [1 - (x-l)/v] - 0 ; 

where fj = number of x f s which exceed j ; and 

a _ 
P - x-1 
v 

13. Shifted Negative Binomial 

p(x) = ^ V+X i 2 ) pV J1 " P,X " 1 ; 

where 0<p<l, v>0 for x = 1,2,... 
The momemt estimates of p and v are given by 

A 

p = x-1 ; and 

s 2 

v = p(x-l) 

1-p 

The maximum likelihood estimator of v is given by the root 

of the equation 

In [1 + (x-l)/v] - ^yr F-; = 0 , 

"7v*j-2) 

where Fj = proportion of x f s # which are greater than 
or equal to j ; and 

A A 

P = -=r-g 

(x+v-l) 

14. Shifted Generalized Negative Binomial 

P(x) = v rrv+flfx-m a**" 1 (i- a )V+£(x-i)-(x-l) 

(x-1)! r[v+/3(x-l)-(x-l)+l] 
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where 0 < a < l, |a/3| < 1, v > 0 for x = 1,2, ... 

The moment esimates of a, /}, and v are given by 

a = l - k& + (h& 2 - l) 3 * ; 

/9 = 1/& [1 - ((5c-1)(1-4)/m 2 )*] ; and 

v = (x-1) (l-4fi)/& 

where A = -2 + [(x-l)^ 3 - 3 u 2 lM. 

(x-1)m 2 3 

and /ii 2 and ^3 are tne second and third central moments 
respectively. 
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